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Abstract. We consider the production of resonant sleptons via R p followed by gauge decays to a charged 
lepton and a neutralino which then decays via R p . This gives a signature of two like-sign charged leptons. 
We find a background at run II of 0.14 ± 0.13 events with an integrated luminosity of 2fb _1 . This enables us 
to probe R p couplings of 2 x 10 for slepton mass of 100 GeV and up to slepton masses of 300 GeV for R p 
couplings of 10 -2 . 

INTRODUCTION 

Sleptons can be produced on resonance at the Tevatron via the fL p LjQjDk term in the superpotential. The 
slepton can then decay again via the K.p-Yukawa coupling. This has been considered elsewhere [1-3]. It can also 
decay via the standard gauge decays of the MSSM [1] which we consider here. There are two possible gauge decays 
of the charged slepton, i.e. 

(1) 

4 - "iX + - (2) 

We focus on the decay of the slepton to a neutralino and charged lepton. The neutralino in turn will then decay via 
the same Jt p -opcrator LjQjDk, 

x° -» {£r + Uj + d k ; vi + dj + d k }. (3) 

Since the neutralino is a Majorana fcrmion it can also decay to the charge conjugate final states, with equal proba- 
bility. In spirit, this is similar to the HERA process considered in [4]. 

The tree-level Feynman diagrams for the slepton production and the neutralino decay are shown in Fig. 1 and Fig. 
2, respectively. Due to the Majorana nature of the neutralino, we have a signature of two like-sign charged leptons. 
In the following we shall consider only electrons or muons, i.e. we focus on the operators L QjDk and LjQjDk- 
We expect these leptons to have high transverse momentum, pt, and be well isolated whereas the leptons from the 
Standard Model backgrounds have lower pt and are also poorly isolated. We therefore hope that this signature can 
be seen above the background if we apply isolation and pr cuts. 



BACKGROUNDS 

In the following we combine the backgrounds for both electrons and muons. The main backgrounds to this like-sign 
dilepton signature are as follows 

1. bb production followed by the production of at least one B° s meson, which undergoes mixing. If the two 
b-quarks in the event decay scmi-lcptonically this gives two like-sign charged leptons. 

2. tt production followed by t — > W + b — ► e + i7 e b, and t — > W~b — ► qqb — > qqW + c — > qqe + i7 e c. 

3. Single top production (s and t channel) followed by semi-lcptonic decays of the top and the B-mcson produced 
after hadronization. 
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FIGURE 1. Production of x°£+. 




FIGURE 2. LQD decay of \ • The neutralino is a Majorana fermion and decays to the charge conjugate final state as well. 




FIGURE 3. Effect of the Parton-Level Cuts. 



4. Non-physics backgrounds from fake leptons and charge misidentification. There are also backgrounds due to 
the production of weak boson pairs, i.e. WZ and ZZ, where at least one of the charged leptons is not detected 
[5]. These require a full simulation including the detector. We do not consider them here. 

We use HERWIG 6.0, [6-8], to simulate these background processes. The program includes the computation of 
the supersymmetric spectrum and the MSSM decay branching ratios from the ISASUSY program [9]. Due to the 
high cross section for the production of bb it was necessary to impose a parton-level cut of 20 GeV on the pt of 
the b and b to enable us to simulate a sufficient number of events. In Fig. 3 we show the distribution of events 
(using the full Monte Carlo simulation) as a function of the (parton-level) pt of the bottom quark for two different 
values of the lepton pt cut. We did not simulate any events for which the pt of the bottom quark was below 20 GeV 
since the cross section is too large. If we extrapolate using Figs. 3a, b to lower b-quark pt we can see that for a 
lepton pt cut of 20 GeV, Fig. 3b, our approximation should be good, i.e. we expect the area under the curve for 
pr(b) < 20 GeV to be negligible. For pt{C) > 15 GeV, Fig. 3a, we would still expect a significant number of events 
at 15 GeV < PT(b) < 20 GeV. Besides the parton-level cut, we forced the B-mesons to decay semi-leptonically. 
This means we neglect the production of leptons from the decay of charmed mesons which should also be a good 
approximation as we expect the leptons produced from these decays to be poorly isolated. 



TABLE 1. Summary of the Background Simulation 



Process 


Cross section 


No. of Events 


Expected No. of Events after cuts, 




before Cuts /nb 


simulated 


for 2fb~ 1 luminosity 


bb mixing 


(9.3 ± 2.3) x 10 4 


8.3 x 10 6 


0.12 ±0.12 


tt 


6.81 ±0.31 


2.0 x 10 5 


0.02 ±0.02 


single top 


1.55 ±0.12 


3.5 x 10 4 


0.00 ± 0.03 


Total 


0.14 ±0.13 



Fig 4a. Effect of an Isolation Cut of 3 GeV Fig 4b. Effect of an Isolation Cut of 5 GeV 
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FIGURE 4. Effect of the lepton isolation cut on the bb background. The upper curve is the full bb background and the 
lower curve is obtained after imposing the isolation cut. 



Figure 4 displays the effect of the lepton isolation cut on the bb background for two different values. The effect of 
the isolation cut on the tt and single top backgrounds is shown in Fig. 5. As can be seen in Figs. 4, 5, by imposing 
an isolation cut of 5 GeV and a cut on the px of the leptons of 20 GeV the background can be almost eliminated. 
Table 1 shows the backgrounds with a pr cut on the leptons of 20 GeV and an isolation cut of 5 GeV. We have used 
the leading-order cross section for the bb and single top backgrounds and the next-to-leading order cross section, 
with next-to-leading-log resummation, from [10] for the tt cross section. In both cases the error on the cross section 
is the effect of varying the scale between half and twice the hard scale, and the error on the number of events is the 
error in the cross section and the statistical error from the simulation added in quadrature. Realistically we cannot 
reduce these statistical errors due to the large number of events we would need to simulate. We have implemented 
the full hadronization using HERWIG 6.0. 

With these cuts and using Poisson statistics, a 5a fluctuation of the total background corresponds to 4 events 
with an integrated luminosity of 2 fb _1 . Hence we consider 4 signal events to be sufficient for a discovery of the new 
R p signal process. 



SIGNAL 

To simulate the signal and the effect of the cuts, we modified HERWIG 6.0, [6-8], to include the production 
process, the MSSM decay of the slepton, and the R p decay of the neutralino. The decay rate of the neutralino 
and its branching ratios were calculated in the code and a matrix element for the neutralino decay [11,12] was 
implemented in the Monte Carlo simulation. 

We use the program to estimate the acceptance of the signal process, i.e. the fraction of the like-sign dilepton 
events which pass the cuts multiplied by the branching ratio to give a like-sign dilepton event. Fig. 6 shows the 
acceptance for two different SUGRA points, with an isolation cut on the leptons of 5 GeV and a cut pr{l) > 20 GeV. 
As can be seen in Fig. 6b, the acceptance drops in two regions. For lower values of M , the slepton is not much 
heavier than the neutralino. The charged lepton from the decay of the slepton is then quite soft and gets rejected 
by the pr cut. For large values of M the slepton is much heavier than the neutralino. The neutralino then gets a 
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Fig 5a. Effect of an Isolation Cut of 5 GeV 
on the top background 
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Fig 5b. Effect of an Isolation Cut of 5 GeV 
on the single top background 
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FIGURE 5. Effect of the Isolation Cut on the tt and single top Backgrounds 
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FIGURE 6. Acceptance from Monte Carlo Simulation for two different SUGRA points. 



significant boost from the slepton decay. The neutralino decay products are folded forward in the direction of this 
boost causing the event to be rejected by the lepton isolation cut. 

To estimate the acceptance properly we need to run a scan of the SUGRA parameter space using the Monte Carlo 
event generator. This still remains to be done. To give some idea of what range of couplings and masses we may 
be able to probe instead we assume an acceptance of 10% using the same cuts as before. We can then estimate the 
range of couplings which may be accessible. 

As can be seen in Figs. 7, 8 the production cross section for A' 211 = 10~ 2 is sufficient to produce a signal which is 
more than 5<r above the background for large regions of the SUGRA parameter space. In some regions where the 
neutralinos become Higgsino-like the cross section drops. The cross section also drops as we approach the region 
where the neutralino is heavier than the slepton and the resonance becomes inaccessible. 

We focused on the coupling X' 211 because the experimental bound on A' m from neutrinoless double beta decay 
is very strict [13,14]. The bound on A' m weakens as the squark mass squared and for squark masses above about 
300 GeV (which we expect in the SUGRA scenario for the heavier slepton masses) A' m w 10~ 2 is experimentally 
allowed and our analysis thus applies for this case as well. \' 2n w 10~ 2 is well within the present experimental 
bounds [14]. 

In these figures we see that we are sensitive to slepton masses up to 300 GeV for couplings of 10~ 2 . The production 
cross section scales with the square of the coupling. For slepton masses around 100 GeV, just above the LEP limits, 
we can thus probe couplings down to about 2 x 10~ 3 . 
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FIGURE 7. Cross section for tan/3 = 2 and A' 211 = 1CF 2 . The long-dashed horizontal line corresponds to the discovery limit 
of 4 signal events for 2fb~ 1 . The signal cross section after cuts is given by the solid and dot-dashed curves for two slepton 
masses. 
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FIGURE 8. Cross section for tan/3 = 10 and A' 211 = 1CT 2 . 



CONCLUSION 

We have performed an analysis of the physics background for like-sign dilepton production at run II and find 
that with an integrated luminosity of 2fb _1 , a cut on the transverse momentum of the leptons of 20GeV and an 
isolation cut of 5 GeV the background is 0.14 ± 0.13 events. This means that 4 signal events would correspond to a 
5<7 discovery, although in a full experimental analysis the non-physics backgrounds must also be considered. 

Using a full Monte Carlo simulation of the signal including a calculation of the neutralino decay rate, its partial 
widths and a matrix element in the simulation of the decay we found that the acceptance for the signal varies but 
for a reasonable range of parameter space is 10% or greater. 

When we then look at the cross section for the production of x°^ + we find that we can probe couplings of 
2 x 10~ 3 for slepton mass of 100 GeV and up to slepton masses of 300 GeV for couplings of 10~ 2 , and higher 
masses if the coupling is larger. 
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